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1 Introduction 

Very shortly after the discovery of asymptotic 
freedom 1 it became apparent that, as a con- 
sequence, at high temperature and/or at high 
density quarks and gluons would also become 
deconhned 2 , leading to a phase transition 
from confined hadronic matter to an uncon- 
fined phase. This was studied in subsequent 
years and since the early 1980ies this phase 
is called the Quark-Gluon Plasma (QGP). 

The conditions for this phase transition 
were studied in lattice QCD and state of the 
art calculations 3 obtain as critical temper- 
ature for the phase transition for two light 
and one heavier quark flavors a value for the 
critical temperature of T c = 173 ± 15 MeV 
and for the critical energy density of e c = 
0.7± 0.2 GeV/fm 3 . It is believed since many 
years that in collisions of heavy atomic nu- 
clei at high energies such conditions should 
be reached. This motivated an experimental 
program starting simultaneously in 1986 at 
the Brookhaven AGS and at the CERN SPS, 
initially with light projectile nuclei such as Si 
and S and from 1992 and 1994, respectively, 
with Au and Pb projectiles. The experi- 
mental results from this program prompted 
a press release from CERN 4 in February 
2000 stating that the combined results from 
the experiments proved that a new state of 
matter other than ordinary hadronic matter 
had been created in these collisions, in which 
quarks were 'liberated to roam freely'. The 
experimental results were clearly not recon- 
cilable with the known hadronic physics and 
it could be estimated that the critical temper- 
ature had been exceeded in the early phase of 



the collision by about 20-30 % and the criti- 
cal energy density by somewhat more than a 
factor 2. On the other hand, from those data 
nothing could be said yet that would charac- 
terize the properties of the new state of mat- 
ter. Hence, at that time the term QGP was 
not used for the new state of matter. 

In the summer of 2000, RHIC as a ded- 
icated collider for heavy ions started oper- 
ation with two large experiments, PHENIX 
and STAR, and two smaller experiments, 
BRAHMS and PHOBOS. In the first 3 years 
of operation data for Au + Au collisions with 
an integrated luminosity of 85/ fib, for p + p 
collisions with 2/pb, and for d + Au collisions 
with 25/nb were collected and a summary of 
the results was recently published in a special 
issue of Nuclear Physics A by all four heavy 
ion experiments 5 ' 6 > 7 ' 8 . In the 2004 Au + 
Au run the 1 /nb level was exceeded and data 
start to appear from this run. Here I will rely 
mostly on published data and review some of 
the key observations from the first 3 years in- 
cluding only a few of the still preliminary first 
run4 observations. 

2 Experimental Results 

2.1 Hadron Production and Statistical 
Models 

Hadron yields have been measured for a large 
range of species at the AGS, SPS and at 
RHIC. It was realized already for many years 
that the data for central collisions of heavy 
nuclei can be rather accurately reproduced 
by calculations for a chemically equilibrated 
system in terms of a grand canonical ensem- 
ble (a review and complete set of references 
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Figure 1. Hadron yield ratios measured at RHIC in comparison to calculations within a statistical model 
based on a grand canonical ensemble (updated version of 10 , taken from 9 ). 



can be found in 9 ). For the lower RHIC en- 
ergy of -\/s = 130 GeV the data are final and 
published and for 200 GeV data are emerg- 
ing currently. Figure 1 shows experimental 
yield ratios from all four RHIC experiments 
in comparison to a statistical model fit. 

In the calculations, there are two free fit 
parameters, the temperature and the baryon 
chemical potential. For top RHIC energy the 
temperature is fitted as 177± 5 MeV, prac- 
tically unchanged from ^/s = 130 and 17.3 
GeV; the baryo-chemical potential is drop- 
ping continuously with increasing beam en- 
ergy reflecting an increasing transparency of 
the nuclei at higher energies and an increas- 
ing dominance of baryon-antibaryon produc- 
tion. This is shown in Figure 2 where results 
of statistical model fits at various beam ener- 
gies are summarized and shown together with 
recent results from lattice QCD 12 . 

It appears that from top SPS energy 
upwards the temperature at which hadro- 
chemical equilibrium is achieved is not chang- 
ing anymore and practically coincides with 
the lattice QCD prediction for the critical 
temperature, while at lower beam energies it 
is falling. At y/s = 8.8 GeV it is only 148 
± 5 MeV. The strangeness suppression that 
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Figure 2. Phase diagram of nuclear matter in the 
temperature - baryon chemical potential plane. Ex- 
perimental points for hadro-chemical freeze-out are 
shown together with a recent lattice QCD calcula- 
tion 12 and a curve of constant total baryon density. 
Figure from 11 . 
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Figure 3. Left: Elliptic flow coefficient V2 as function of pt for different particle species . Together with 
the experimental data results from a hydrodynamics calculation including a phase transition are shown 30 . 
Right: Experimental data for pions and protons at a lower RHIC energy 26 . Also shown are hydrodynamics 
calculations 30 with and without phase transition from QGP to hadronic phase. Figure from 7 . 



is well established for pp and e + e~ collisions 
appears to be completely lifted. This leads 
to an enhancement in the yields of particu- 
larly multistrange hadrons in heavy ion col- 
lisions as compared to pp results. For the 
Omega baryon at SPS energy this enhance- 
ment is 13 a factor 17. How hadrons like the 
can be equilibrated on the time scales of the 
nuclear collision has been a puzzle for sev- 
eral years and there is consensus that with 
two-body collisions and the known hadronic 
cross sections this is not possible 14 > 15 > 16 . A 
possible explanation has been presented re- 
cently 14 . In the direct vicinity of the phase 
transition the densities of particles are rising 
very rapidly due to the increase of degrees of 
freedom by more than a factor of 3 between a 
hadron gas and a QGP. At these high densi- 
ties multi-hadron collisions become dominant 
and can drive even the f2 yield into equilib- 
rium in a fraction of a fm/c. Conversely, 
already 5 MeV below the critical tempera- 
ture the densities are so low that the system 
falls out of equilibrium and the yields can- 
not follow anymore a decreasing temperature. 
Therefore the authors of 14 conclude that the 



rapid equilibration is a direct consequence of 
the phase transition from QGP to hadronic 
matter and that, at least at high beam ener- 
gies, the chemical equilibration temperature 
is a direct experimental measure of the criti- 
cal temperature. 

2.2 Elliptic Flow 

Momentum distributions in three dimensions 
are analyzed with transverse coordinates rel- 
ative to the reaction plane of the collision 
spanned by the impact parameter vector and 
the beam direction and a decomposition in 
terms of Fourrier coefficients is performed. 
Already at the Bevalac sizeable anisotropies 
were observed for heavy colliding nuclei. In 
particular, the quadrupole coefficient V2 was 
found to be negative, explained by shadow- 
ing of the emitted particles by the target and 
projectile spectator remnants 17 . At AGS en- 
ergies a sign change was observed 18 by E877, 
i.e. the momentum spectra were harder in 
the reaction plane than perpendicular to it. 
The interpretation used a prediction from hy- 
drodynamics 19 that, for semiperipheral col- 
lisions, in the early phase of the collision the 



For Publisher's use. 




Pj (GeV/c) pj (GeV/c) 

Figure 4. Neutral pion transverse momentum spectra measured by PHENIX in peripheral (left) and central 
(right) Au + Au collisions (stars) together with pp data from the same experiment scaled with the number 
of binary collisions (circles) 37 > 38 . Yellow band: Normalization uncertainties of the pp data. Black line: NLO 
pQCD calculation. Figure taken from 39 . 



pressure gradient was larger in this direction 
due to the excentricity of the nuclear overlap 
region . From the hydrodynamic evolution 
it would follow, that this anisotropy in pres- 
sure gradient would evolve with time into an 
anisotropy in momentum space, driven by the 
initial condition and the equation of state of 
the expanding system. This was confirmed 
by a microscopic analysis within a transport 
model 21 . From this the name 'elliptic flow' 
originated for the quadrupole coefficient v 2 . 

At the higher SPS energy growing posi- 
tive coefficients V2 were found 22 ' 23 > 24 and the 
sign change was traced to occur 25 at beam 
momenta per nucleon of about 4 GeV/c. 
At RH1C energies very large values of v 2 
were observed 26 . 27 . 28 > 29 j typically about 50 % 
above SPS top energy results. 

This was studied differentially for differ- 
ent hadronic species and as function of p t as 
shown for data from STAR in Figure 3. ft is 
observed that for more massive hadrons the 
rise of v 2 starts at larger values of p t . For 
the first time there was quantitative agree- 

a The use of hydrodynamics to describe the dynamics 
of a hadronic collision goes back to the 1950ies 20 . 



mcnt with hydrodynamic calculations 30 > 31 
in terms of pt and hadronic species depen- 
dence, as also shown in Fig. 3. These hy- 
drodynamic calculations also reproduce the 
overall features of the pt spectra of differenct 
hadrons, although in details there are devia- 
tions stemming from the different treatment 
of the hadronic phase and freeze-out (see Fig. 
20 in 8 and references there). It is common 
to all the hydrodynamics calculations that, 
in order to reproduce the data, a rapid ini- 
tial equilibration on a time scale faster than 
1 fm/c is required 30 < 31 < 32 . 

At pt above 2-3 GeV/c, where hydrody- 
namics should no longer hold as a theoretical 
description, another type of scaling was dis- 
covered 33,7 : dividing both v 2 and pt by the 
number of constituent quarks in a hadron all 
results match rather well even including mul- 
tistrange baryons. It was realized that an 
old idea of quark coalescence 34 could be the 
underlying physics 35 and indeed calculations 
based on the assumption of coalescence of va- 
lence quarks during hadronization of a QGP 
reproduce this feature rather well 36 . 
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2.3 High Momentum Suppression 

One of the highlights of the RHIC experimen- 
tal program is the observation of a strong 
suppression in the production of hadrons 
at high transverse momentum when com- 
pared to pp collisions. Figure 4 shows the 
Pt spectrum of neutral pions in Au + Au col- 
lisions as compared to a measurement in pp 
in the same experiment and at the same en- 
ergy 37 > 38 . The pp spectrum compares well 
with a calculation in NLO pQCD. In order 
to compare, the pp spectrum has been scaled 
with the number of binary nucleon-nucleon 
collisions in a Au + Au collisions at a given 
centrality. The number N co ii of binary col- 
lisions is given by the collision geometry - 
measured in the data with some resolution 
-, the well known nuclear density distribu- 
tion, and the inelastic pp cross section. The 
collision centrality in Au + Au collisions is 
characterized by the fraction of the geometric 
cross section for which events have been se- 
lected, which is related to the impact param- 
eter. The yellow bands in Fig. 4 reflect the 
systematic uncertainty in this scaling. One 
can observe that for peripheral collisions pp 
and Au + Au collisions agree very well, while 
in central collisions the Au + Au spectrum is 
significantly suppressed. 

This is better visualized by building the 
ratio Raa between the Au + Au pt spectrum 
and the pp spectrum scaled with N co ii as 
shown in Figure 5. All four RHIC experi- 
ments observe a suppression by about a fac- 
tor of five for p t larger than 4 GeV/c. Since 
not all experiments measure neutral pions, 
the ratio is shown here for charged hadrons, 
but at large pt the data for all hadron species 
merge. At low p t the ratio Raa is expected 
to be below one because there, due to the 
dominance of soft processes, the appropriate 
scaling is with the number of participants, i.e. 
nucleons in the nuclear overlap region. It is 
expected that this ratio should rise as hard 
scattering becomes dominant and, in fact, 
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Figure 5. Ratio Raa of the pt spectrum for cen- 
tral Au + Au collisions normalized to the pp spec- 
trum scaled with the number of binary collisions for 
charged particles from all four RHIC experiments. 
Figure from 40 . 



due to the well known Cronin enhancement, 
in the region of 2-6 GeV/c values above one 
are expected. Contrary to this expectation 
the data show a suppression. 

The suppression is not unexpected. It 
was predicted that in a medium with high 
parton density the radiative energy loss of 
a quark or gluon should be strongly en- 
hanced 41 > 42 J leading to a very effective ther- 
malization of jets in a hot color charged 
medium. Calculations employing a large ini- 
tial gluon rapidity density of about 1100 can 
account 43 for the data at top RHIC energy. 
The beam energy dependence of the Raa ra- 
tio was presented recently by d'Enterria 44 
and it appears that the suppression evolves 
in a very smooth way from top SPS energy 
onwards. The RAA(pt = 4GeV/c) values are 
shown in Figure 6 for the top SPS energy 
and three RHIC energies. Already the val- 
ues of about 1.0 measured at the SPS repre- 
sent a slight suppression as compared to the 
normal Cronin enhancement 44 . Going from 
y/s^ = 17.3 to 62.4 to 200 GeV the gluon ra- 
pidity density needed to reproduce the data 
grows 43 from 400 to 650 to 1100. An al- 
ternative formulation of this in medium sup- 
pression is by increasing and large opacities 
of the medium traversed 45 . 



For Publisher's use 



> 

u 
(5 



_ Suppression at p T = 4 GeV/c: 

■ Pb+Pb -> n°+X 0-7% central [WA98] 
A Pb+Au -* T^+X 0-5% central [CERES] 

S+Au -» 7i°+X 0-8% central [WA80] 
• Au+Au -» n°+X 0-1 0% central [PHENIX] 



0.5 
0.4 

0.3 
0.2 




X.N.Wang jet quenching: 

Non-Abelian energy loss: A E /A E q = 9 
"Non-QCD" energy loss: A E = A E 



0.1 r 



PHENIX Au+Au {central collisions): 
■ Direct y 
A af Preliminary 



20 40 60 80 



100 120 140 160 180 200 220 

\ls»,(GeV) 



Figure 6. Suppression factor Raa f° r pions as func- 
tion of beam energy at a fixed pt value of 4 GeV/c. 
The yellow line shows a calculation 45 including par- 
ton energy loss in a medium with high parton density. 
Figure from 44 . 



The proof that this is really a final state 
effect probing the properties of the medium 
traversed by the parton is given by the ob- 
servation that in d + Au collisions in the 
same experiments no suppression is seen, 
but rather the expected Cronin enhance- 
ment 46,47,48,49_ 

Direct photons were measured by 
PHENIX in pp and Au + Au collisions at top 
RHIC energy 50 ' 51 . The pp spectra are rather 
close to a NLO pQCD calculation 52 . The Au 
+ Au photon spectra are within errors con- 
sistent with the scaled pp result and hence 
the expectation from NLO pQCD. For all 
centralitics they do not show any significant 
suppression. This is shown in Figure 7 where 
also the most recent neutral pion results from 
run4 53 extending out to p t = 20 GeV/c are 
displayed. It is remarkable that the suppres- 
sion of the pion p t spectrum remains practi- 
cally constant over a large range in p t from 
4 to 20 GeV/c, close to the predicted behav- 
ior for a medium with initial gluon rapidity 
density of 1100 (see Fig. 7). 

The high initial gluon densities corre- 
spond to an initial temperature of about 
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Figure 7. Preliminary PHENIX results for the sup- 
pression factor Raa ou t to high pt for n° and rj 
mesons together with a calculation based on a high 
gluon rapidity density 43 . Also shown are the results 
for direct photons. Figure from 53 . 



twice the critical temperature and to initial 
energy densities eo = 14 - 20 GeV/fm 3 well 
in line with the initial conditions needed for 
the hydrodynamics calculations to describe 
spectra and elliptic flow (see previous section) 
and bracketed by the estimates based on the 
Bjorken formula and transverse energy pro- 
duction. 

The observed high p t suppression pattern 
is different for different hadronic species 5 ' 8 ' 7 . 
In particular, a pattern appears where at in- 
termediate values of p t of 2-6 GeV the sup- 
pression of baryons is significantly weaker 
than that of mesons. The proton/pion or 
also the A/K° ratios peak at values 1.5-1.6 
for p t = 3-4 GeV/c, close to the ratio 3/2 ex- 
pected in quark coalescence models. 

Parton thermalization is displayed in a 
very clean way be recent results of the STAR 
collaboration 54 . Evaluating the mean trans- 
verse momentum in a cone opposite to a high 
Pt trigger particle as a funtion of centrality, a 
gradual decrease for more central Au + Au 
collisions is observed and in the most central 
collisions a value very close to the inclusive 
mean p t is reached (see Figure 8). 

In azimuthal correlations of two high 
Pt particles it was seen that the away-side 
peak disappears in central Au + Au collisions 
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Figure 8. Mean transverse momentum in a cone of 
an opening angle of one radian opposite to a leading 
particle as a function of the collision centrality; dis- 
tributions are for two ranges of leading particle trans- 
verse momentum as well as for the inclusive particle 
distribution. Data and figure from STAR 54 . 
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for a choice of trigger p t of 4-6 GeV and pt of 
the correlated particle of 2-4 GeV/c 55 . In 
pp, d + Au and peripheral Au + Au collisions 
a clear peak opposite to the trigger particle 
is observed in the same type of correlation, 
also measured by STAR 55 ' 48 . Recently, it 
was shown that the effect is very strong in 
case the away-side jet is emitted out of the 
reaction plane and much weaker for emission 
in the reaction plane 56 as displayed in Fig- 
ure 9. This supports the strong correlation 
of the suppression with the length of matter 
traversed by the parton. 

When lowering the pt cut on the corre- 
lated hadron, a very broad structure appears 
on the side opposite to the trigger parti- 
cle. This was shown by STAR 54 for a cut 
on the correlated hadron of p t = 0.15 - 4 
GeV/c. This calls to mind a similar obser- 
vation at SPS energy by CERES 57 where for 
a condition p t > 1.2 GeV/c for both particles 
also very strong broadening of the away-side 
structure with increasing collision centrality 
in Pb + Au collisions was observed. Recent 
data 58 from PHENIX display a tantalizing 
feature as shown in Figure 10: For a trigger 
particle pt of 4-6 GeV/c and a correlated par- 
ticle pt of 1.0 - 2.5 GeV/c the away-side peak 
seen in peripheral Au + Au collisions devel- 
ops actually into a hole at A(f> = it for more 
central collisions while a very broad peak ap- 
pears with a maximum at Atp = tt — 1 as 
can be seen in Fig. 10. A suggestion has 
been made that this could be the Mach cone 
due to the sonic boom of the quenched jet. 
A parton traversing a quark-gluon plasma 
with velocity larger than the velocity of sound 
in the QGP (a/1/3 for an ideal gas) would 
radiate only up to a cone angle of about 



Figure 9. Azimuthal angle correlation between a high 
pt trigger particle (4-6 GeV/c) with all particles in 
the window pt = 2-4 GeV/c for intermediate central- 
ity Au + Au collisions. The distribution is shown for 
both particles in a ± 45 degree window around the re- 
action plane orientation (in-plane) or a same window 
perpendicular to it (out-of-plane). Figure from 56 . 



would have far reaching consequences since it 
would be an observable linked directly to the 
speed of sounds of the quark-gluon plasma 
and thereby its equation of state. It remains 
an experimental challenge to establish an ac- 
tual cone topology in two dimensions. 
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Figure 10. Azimuthal correlations of a leading par- 
ticle of pt = 4-6 GeV / c and any particle with pt = 
1-2.5 GeV/c for different centralities of a Au + Au 
collision. Data and figure from PHENIX 58 . 

2.4 Charm Quarks and Quarkonia 

Open charm has been measured indirectly 
from the inclusive electron p t spectra after 
subtracting known contributions from pho- 
ton conversions and light hadron decays by 
PHENIX 61 . The spectrum remaining after 
subtraction is dominated b by open charm 
and beauty contributions. Recently results 
for an elliptic flow analysis were shown 62 
of the electrons dominantly from open charm 
decays. There is a significant nonzero value 
in the pt range 0.4 - 1.6 GeV/c. This is con- 
firmed by preliminary STAR data 63 that ex- 
tend the overall transverse momentum cover- 
age by adding the range p t = 1.5 - 3.0 GeV/c. 
Together, the data paint a consistent picture 
that indeed the electrons from open charm 
decay exhibit elliptic flow, i.e. follow the 
collective motion of the light quarks. This 
would imply that the charm quark thermal- 
izes to a significant degree. Note that this is 

b A possible contribution to the electron spectrum 
from the Drell-Yan process cannot be ruled out at 
present, though. 
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Figure 11. Raa suppression factor for electrons dom- 
inantly from open charm and open beauty decay for 
Au + Au collisions at top RHIC energy for different 
collision centrality (see text) measured by PHENIX. 
Figure from 64 . 

a necesarry prerequisite for any formation of 
charmed hadrons by statistical hadronization 
(see below). 

In that case also jet quenching should be 
observed for charmed hadrons. Indeed, in 
still preliminary data it was shown recently 
that electron spectra, after the subtraction 
of contributions from conversion and light 
hadron decays, show high p t suppression for 
central Au + Au collisions 64 . The Raa fac- 
tor drops practically as low as for pions at 
Pt of 4 GeV/c, i.e. to values of about 0.2. 
In a recent publication 65 the suppression for 
electrons from D meson decay was studied 
for different transport coefficients. The pre- 
liminary RHIC data would be consistent with 
a calculation using a transport coefficient of 
14 GeV 2 /fm (see Fig. 2 of 65 ), at the upper 
end of the range needed to reproduce Raa for 
pions. This is very surprizing, in particular 
also in view of the fact, that at pt of about 
4 GeV/c also the contribution of b-quarks 
to the electron spectrum should become size- 
able. 

At top SPS energy, for central Pb + Pb 
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Figure 12. J/ip yield in Au + Au and Cu + Cu colli- 
sions at top RHIC energy normalized to the measured 
result for pp collisions, scaled with the number of in- 
elastic collisions. Results are shown for decays into 
electron and muon pairs at mid- and forward rapidi- 
ties, respectively. Figure from 53 . 



collisions a significant suppression of J /ip pro- 
duction was observed in the NA50 experi- 
ment 66 > 67 . This suppression is compared to 
the socalled 'normal' nuclear absorption seen 
also in pA collisions. From analysis of all 
pA data, a cross section for normal nuclear 
absorption of 4.1 ±0.4 mb was extracted 67 . 
To this normal nuclear absorption all results 
from heavy ion collisions can be compared. It 
turns out that S + U data as well as data from 
peripheral Pb + Pb collisions agree with this 
normal nuclear absorption curve. For trans- 
verse energies above 40 GeV or a length of 
nuclear matter seen by the 3/tp of L > 7 fm 
the points from Pb + Pb collisions fall in- 
creasingly below this normal nuclear absorp- 
tion curve. Theoretically, the suppression can 
be explained by disappearance of the 3/tp (or 
possibly only the charmonia states that feed 
it) in a hot colored medium or by interaction 
with comovers (chiefly pions), albeit with a 
very large density of more than 1/fm 3 , i.e. 
a value not deemed achievable for a hadron 
gas. 

The first results for 3/tp production in 



central Au + Au collisions at RHIC energy 
came very recently from PHENIX 53 (run4); 
they are displayed in Figure 12. As compared 
to pp as well as d + Au collisions there is a 
significant suppression. The suppression is, 
however, rather similar to the one observed 
at SPS. This is in contrast to some predic- 
tions 68 that in central Pb + Pb collisions 
at RHIC only 3/tp mesons from the corona 
should survice, i.e. order of 5 % of the normal 
unsuppressed yield. The actually observed 
yield by far exceeds this expectation. This 
could be seen as indication, that at RHIC the 
energy density in the QGP is not yet high 
enough to dissolve the 3/tp but rather only 
enough to dissolve higher cc states. Recent 
results from high temperature lattice QCD 
indicate that the 3 /tp bound state may only 
disappear 69 in the vicinity of T = 2 T c . 

A maybe more interesting alternative 
has been proposed 70 : Even if the initially 
formed cc pairs are completely dissociated 
in the hot QGP, at hadronization charmed 
hadrons may form in a statistical fashion 
by the same mechanism described above for 
hadrons involving up, down, and strange va- 
lence quarks. This includes also the forma- 
tion of charmonia and it was pointed out in 70 
that the formation probability oi 3/tp mesons 
would grow quadratically with the cc rapid- 
ity density. Such increased reformation of 
3/tp by statistical hadronization could possi- 
bly account for a suppression apparently not 
much stronger than at SPS. 

Figure 13 shows the prediction of 71 to- 
gether with the recent and still preliminary 
RHIC data for 3/tp — > e + e~. Currently the 
main uncertainty in this consideration is the 
overall charm production yield, which enters 
quadratically. A calculation is shown using 
a cross section per unit rapidity of 80 /A> 
(400 /ib integrated) as expected from NLO 
pQCD 68 . The overall charm production 
cross section at RHIC energy has so far been 
measured indirectly by PHENIX 61 for 
= 130 and 200 GeV from the inclusive elec- 
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Figure 13. J/?/' yield per participating nucleoli for Au 
+ Au collisions at top RHIC energy compared to the 
yield expected from statistical recombination. Data 
from 53 , statistical hadronization prediction from 71 . 
For statistical hadronization a standard interval of A 
y = 1 is used. Calculations are shown for a cc cross 
section per unit rapidity of 80 fib, of 160 fib. Also 
shown: a calculation for 80 fib and A y = 2.0. 



tron spectra in the way described above for 
Au + Au collisions and, at full RHIC energy, 
also for d + Au and pp collisions 72 . It is 
found that the integrated charm cross sec- 
tion, when scaled with the number of binary 
collisions, agrees for all three collision sys- 
tems. The value is about 50% above the NLO 
pQCD calculation 68 but agrees within er- 
rors. On the other hand, in STAR, D mesons 
have been reconstructed via their hadronic 
decay to Kir in d + Au collisions and a charm 
cross section per nucleon nucleon collision has 
been extracted 73 . It is twice as large as the 
PHENIX value by nearly two standard devi- 
ations. The experimental situation concern- 
ing open charm production needs to be im- 
proved before the J/ip puzzle can be better 
addressed. Only measurements at LHC will 
unambiguously clarify the role of statistical 
hadronisation of charm, since with this mech- 
anism a significant 3/ip enhancement in Pb + 
Pb collisions at LHC was predicted 71 instead 



of suppression. Such an enhancement would 
be an unambiguous signal of deconfinement. 

3 Summary and Outlook 

Hadron yields are found to be in chemical 
equilibrium. For top SPS energy and up this 
can be achieved by multi-particle collisions 
in the direct vicinity of T c and hence the ob- 
served chemical equilibration temperature is 
an experimental measure of the critical tem- 
perature for the phase transition. 

At RHIC energies, spectra and azimuthal 
correlations are quantitatively described by 
hydrodynamics. This requires rapid local 
thermalization and high initial energy den- 
sities more than tenfold above the calculated 
critical energy density for the phase transi- 
tion between hadronic matter and QGP. 

High p t hadrons are suppressed in central 
Au + Au collisions and this is a medium ef- 
fect. Jet quenching in a hot color charged 
medium was predicted, modelling of the data 
with high parton density is successful. There 
are some indications of valence quark coales- 
cence in hadronic observables. 

The observations that lead to the CERN 
press release are confirmed by the RHIC ex- 
periments. Beyond this additional features 
are observed that start to probe the prop- 
erties of the new state of matter. Much 
progress in this direction is expected from the 
high luminosity RHIC data just starting to 
appear and, from 2007, from the LHC heavy 
ion program. 
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